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Abstract 



We discuss the non-relativistic multichannel quark model and describe the 
techniques developed to solve the resulting equations. We then investigate 
some simple solutions to demonstrate how the model unifies meson-meson 
scattering with meson spectroscopy, thereby greatly extending the domain of 
applicability of the naive quark model. In the limits of narrow resonance 
widths and no quark exchange, it reproduces the standard quark model spec- 
troscopy and Breit-Wigner phase description. Outside those limits s-channel 
resonance masses are lowered by their two-meson couplings, the line-shapes 
of wide resonances are significantly altered, and the equivalent Breit-Wigner 
masses and widths show an energy dependence. Because meson-meson in- 
teractions are due to coherent s-channel resonance production and t-channel 
quark exchange (though other interactions can readily be added), the multi- 
channel equations model experimental resonance production and decay in a 
way that the usual eigenvalue equations cannot. 
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I. INTRODUCTION 



The quark model attempts to represent the main features of QCD, the theory of quarks 
and gluons, and has been successfully used to calculate hadron spectroscopy and decays [0. 
However, after decades of study, and numerous extensions, a multitude of discrepancies still 
remain between theory and experiment, particularly in the light meson sector. 

In this paper we construct and solve a multichannel quark model, which we apply to 
meson-meson scattering. We consider only low energy scattering processes, which are res- 
onably soft for all mesons except the vr at the high end of our energy range, so we use 
a non-relativistic formulation. We motivate the model phenomenologically by considering 
strong decay processes, and add terms representing these processes to the naive quark model. 
Numerical solutions, which are described in the Appendicies, yield the S-matrix parameters 
for the free two-meson states and the masses, widths, and two-meson decay predictions for 
the s-channel meson resonances. The standard quark model spectroscopy and the Breit- 
Wigner phase shift are reproduced, in the limits of narrow resonance width and vanishing 
quark exchange. Outside those limits resonance properties are found to depend significantly 
on two-meson couplings, which leads to mass shifts, energy dependent Breit-Wigner param- 
eters, and two-meson bound states. 

The paper is organized as follows: We first review spectroscopy in the naive quark model. 
Next we develop the multichannel equation for the simple case of the elastic scattering of two 
stable mesons, A and B, through an intermediate resonance R. Some qualitative behaviours 
of the solutions are explored, new effects are highlighted, and likely avenues of further inves- 
tigations are noted. We then review the quark exchange process and incorporate it into the 
multichannel model. Finally we apply the model to inelastic multichannel scattering. 

This model has been discussed previously in the square well approximation ||^,^ and, in 
a less complete form in . See also . 

The fortran source code that has been developed to find the numerical solutions to the 
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multichannel equations, and rudimentary intructions for running it, are available over the 
World Wide Web §. 



II. THE NAIVE QUARK MODEL 

One method of predicting the meson spectrum from the naive quark model is to solve 
the eigenvalue problem 



(1) 



for each qq meson resonance R. Here the qq Hamiltonian i^H^_(rgq) defines all the properties 
of R. In the naive quark model we use a phenomenological qq central potential Vn -(r^g), so 



that equation (1) becomes 
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where \4 -('"gg) takes the simple form 
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The eight components of Fg (Fg) are defined by [Fg]° = ^X"^ ([F,]" = -|(A")*), for a = 
1, 8, where the A" are the eight 3x3 SU(3) color matricies. The and Sg are defined by 
[Sg]* = [Sg]* = |o"*, for i = 1,3, where the a are the three 2x2 SU(2) spin matricies. The 
expectation value of the Fg ■ Fg factor reduces to a constant factor of —4/3 for mesons, while 
(Sg ■ Sg) = — 3/4 for spin zero qq systems and +1/4 for spin 1 systems [|l|. 

The linear confinement term (&rgq) is predicted by the flux-tube model and by lattice 
QCD. The adrqq term is a short-range Coulomb piece that dominates the interaction when 
Tgg -C 1 fm. It can also be represented by a short range gaussian potential, which approxi- 



mates the mid- and long-range behaviour of the 1/Tgg interaction, but is finite at rgg = 0. 



C is a constant background potential. The last term is the short-range color magnetic 
hyperfine interaction. It increases (decreases) the energy of systems whose color magnetic 
moments are aligned (anti-aligned). We ignore spin-orbit coupling terms, which normally 
lead to mass splitting of order 10 MeV between, for example, the "^Pq, ^Pi, and ^P2 states. 
In Section |TTl| we will show, for the iqq = 1, Sqq = 1, Kq, Ki, and K2 states specifically, that 
different couplings and different external two-meson states shift the underlying qq resonance 
masses by 10s of MeV. 

If we use separation of variables to write the qq wavefunction as 

V'fi,, (r,,-) = u^^_ (rqq) {9, 0) (4) 

then Uji^-^{rqq) obeys the radial equation 

Hr,, (rqq) (rqq) = {t qq) (5) 

where Hii^^{rqq) is the radial part of the qq Hamiltonian given by 

_(rqq) = -J—-^ + + 1) + ^ + ^ + Vn - (rqq). (6) 

"■qq\ 99/ r) _^ -S o,, _ r ^ ^ "-qq \ 99/ V / 

r'qq ^' 99 ^ H'qq ' qq 

In ^qq is the reduced qq mass, given by fiqq = 2mqm^/{mq + rriq). The second term is the 

centrifugal barrier term, which comes from the three dimensional kinetic energy operator. 

Since we limit ourselves to this one-dimensional equation, we choose to view the barrier 

term as part of the potential, rather than part of the kinetic energy. 

Solutions of equation (|^) yield meson masses and meson wavefunctions ^^^-(rqg) which, 

here, depend on quark flavor, qq separation r^g, spin Sqq, and orbital angular momentum Iqq. 

The light quark sector of this Hamiltonian has ten parameters, m^, md, m^, b, ac, a^, Cg, 

C, ah, and cxh, which are fixed so the model reproduces a small subset of the available data. 

Predictions for the masses of all other light mesons follow. 

As we are using constituent, rather than current quark masses, we impose the constraint 
= rrid, reducing the number of parameters to nine. We also choose ag = 0, so the value 

of Cg becomes irrelevant, and we have seven free parameters. 



In this paper we use m„ = = 0.375 GeV, = 0.650 GeV, b = 0.782 GeV/fm, 
= 0.857, C = -0.581 GeV, ah = 0.840, and df, = 0.700 GeV, which leads to the meson 
masses quoted below. 

In Section fT^ we will add two more parameters that couple qq resonances and two-meson 
states, and the scattering properties of two-meson systems can then be determined. 

The wavefunctions ^^^-(rgg) solving equation (1) are found numerically by finite difference 
methods using the Noumerov technique The details are provided in Appendix 0. In 
the large r region, these numerical solutions are, generally, linear combinations of an expo- 
nentially decaying and an exponentially growing term. The relative amplitude of each term 
depends on the trial energy Ej^. The physical solutions are defined by the additional con- 
straint that Uji{rqq) ^ as — * oo, which requires that the amplitude of the exponentially 
growing factor vanish, and defines the eigenenergies m^j. Details are given in Appendix ^. 

In FIG. 1. we plot the Iqq = 1, Sqq = 1, -ft'j=o,i,2 radial wavefunctions at trial energies of 
El = 1.4299235012 GeV and E2 = 1.4299235014 GeV. (We have no spin-orbit term in the 
Hamiltonian, so the Kfs are degenerate at this point.) The lower energy solution goes to 
+00 as Tqq — ^ 00 while the other solution, which is only 0.2 eV greater in energy, goes to —00 
as Tqq — > 00! The eigenenergy lies somewhere in between these two values, so we can clearly 
use this proceedure to determine the eigenenergies of @ to an irrelevantly high accuracy. 

As a verification of this technique we tested the iterating program on the toy problem of 
two quarks confined by an harmonic oscillator potential. We found that the trial energies 
differ from the exact analytic eigenenergies by less than one part in 10^^, and the numerical 
wavefunctions differ from the exact gaussian solutions by less than one part in 10^ out to 
distances of 8 fm! The harmonic oscillator is a particularly simple potential however, and we 
must always check the solutions to make sure they are independent of the various iteration 
parameters introduced to numerically solve equation (^). 

In FIG. 2. we show the wavefunction of the first radial excitation of the Kj at energies of 
El = 2.047471763 GeV and E2 = 2.047471768 GeV. This energy resolution of 5 eV is coarser 
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than above simply to make the wavefunctions diverge a distances of less than 6.0 fm. (As 
the trial energy approaches the eigenenergy the point of divergence of the trial wavefunction 
moves further from the origin.) The second radial excitation is at 2.5651608366 GeV, and 
the higher radial excitations are easily found this way. 

With the Hamiltonian parameters above, this approach yields meson masses (in GeV) of 
= 0.133, TTiK = 0.494, = 0.775, tuk* = 0.913, = 1.019, = 1.230, rriK, = 1.389, 
= 1.512, rriar = 1.280, rriK, = 1.430, and mf = 1.550, with J=l, 2, and 3, in good 

1 •'J 

agreement with experiment. 

III. S-CHANNEL MESON RESONANCE PRODUCTION 

When solving any qq eigenvalue problem one implicitly assumes that resulting mesons 
have a unique energy, and are therefore stable. However, most hadrons are hadronically 
unstable, with typical widths of 200 MeV. The Heisenberg uncertainty principal gives them 
a lifetime of 

r^J^^ lO'^Hm ^ 3 X lO'^Sec , (7) 
AE ' ^ ^ 

or as long (short) as it takes light to cross a meson. The zero width approximation intrinsic to 

the eigenvalue formulation is, therefore, unrealistic for all the hadronically decaying hadrons. 

When we try to extract the spectroscopy of wide and interfering hadrons, glueballs, and 
hybrids from experimental data, as a means of gathering information about the underlying 
quark-quark interactions, and to test theories of QCD, we find the predictions of equation (^) 
inconsistent with data. There are results, such as the presence of two p' states, which have 
remained unexplained after years of intense research and many extensions of (§). 

The multichannel model we propose scatters two-meson states into two-meson states 
through s-channel resonance production and t-channel quark exchange. The relative wave- 
function of two-meson system, and the wavefunction of s-channel resonance, are found. 
From these the S'-matrix and the resonance parameters are deduced. 
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Consider s-channel meson resonance formation. As an idealized gedanken scattering 
experiment one might think of two mesons, A and B, coalescing via qq annihilation to form 
a single meson resonance R. (We will include the t-channel quark exchange effects below.) 
Some time later R decays, via qq creation, into a final two-meson system CD, where CD 
may or may not represent the same mesons as AB. This is depicted in FIG. 3. 

For the moment we restrict our discussion to the elastic scattering process AB ^ i? — > 
AB. Consider the initial state AB in FIG. 3. The Schrodinger equation, which determines 
the relative AB wavefunction MAs(r4s), is 

Hab {^ab ) Uab {^ab ) = Eab UabItab) , (8) 

where the radial part of the two-meson Hamiltonian HABi^^AB) is given by 



HAB{rAB) = ( T— ^ ~ ^Ab{Lb + 1) I ^ VAB^AB{rAB) + ^Ia + rUg ■ (9) 

^ Mas V O-Vab ^ab 



We have included the meson masses, tua and ntg, in the Hamiltonian so that the energy 
Eab represents the total energy of the AB system in its center of momentum frame. The 
potential Vab^ab will represent other possible non-resonant AB interactions, such as quark 



exchange shown in FIG. 4., and will be discussed in Section 0. For the moment we take 

Vab^ab — 0. 

Equation (||) says that the AB energy operator, Hab, acting on the AB probability 
amplitude, MAs(r4s), equals the availiable energy Eab times UAfl(?"As)- When one looks at 
the transition AB to R in FIG. 3., however, one sees that R removes energy from AB. 
This suggests modifying equation (|D by removing a separation dependent energy factor, 
Ias^h('as) Uqq{rqq), from the right hand side of (H), leading to 

HabUab{^Ab) =EabUab{^Ab) ~ VAB^ai^AB) Uqq^Tqq) . (10) 

A standard model for qq annihilation couples the qq state to the vacuum, which requires 



that the annihilating qq pair be in an igg = 1, Sqq = 1, ^Pq color singlet state |[TO| , p!T|| . It 
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is, however, possible that some annihilation verticies are dominated by the process qq g 
(where this g is a gluon). This would require that the annihilating qq pair have the ^5*1 
quantum numbers of the gluon, constraining them to an igg = 0, Sgg = 1, color octet state, 
and making R a hybrid (qqg) meson. The model we are exploring here can be adapted to 
quantitatively study this fundamental question, though we do not do so here. 

In this paper we choose, as an ansatz for the creation and annihilation verticies, the 
transition potential 
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where g is an unknown universal constant (with dimensions of energy) to be determined 
from some subset of experimental data. C^l^^jj and C/g^^ are, respectively, the spin and 
flavor Clebsch-Gordon factors relating the overlap of the AB and R spin and flavor states, 
and depend on the structure of the qq creation and annihilation verticies. We generally take 
^AB^R — ^AB^R — 1 here, an approximation that needs to be removed before quantitative 
predictions can be realized. Iab indexes the angular momentum of AB. The second unknown 
parameter in ( pH]) is Tq, which sets the range of the interaction. It must be determined from 
flts to experiment but is expected, from naive arguments, to be approximately 0.5 fm. In 
this paper we use Tq = 3.0 GeV^^ ~ 0.6 fm. Ni^^ is a normalization constant deflned by the 
constraint 
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e-i-AB/ro? ^ ^ 1 _ (12) 



Other forms of the transistion potential VAB^R{rAB) are obviously likely [^, and the 



exploration of these possibilities is an important application of the multichannel model. 

Notice that both r^B and Vgg appear in equation (0), so we must relate them. An 
examination of the qqqq geometry (FIG. 5.) an instant before the transition into the qq state 
shows, in the SU(3) limit where m„ = = m^, that = Vgg = 2rAs. This relationship is 
easy to implement in (|10|). (See Appendix p.) 



The qq resonance R in FIG. 3. was described by equation (^). From FIG. 3., however, 
we see that it must be modified to allow for the transfer of the energy from R into the AB 
system, leading to 

By time reversal invariance we require that Vab^r = yR_,AB = Vab^r- 



IV. THE TWO CHANNEL MODEL 



Equations ([T0|) and (|T3]) can be combined into matrix form: 



Hab{^ab^ 



AB^RV AB ) 



AB^r{^Ab) Hji^-{rqq) 



Uab{Xab') 

UR,Arqq) 



(14) 



Uab{Xab^ 

where E^g is the energy available in the AB center of momentum system and takes on values 
greater than + nig Equation is the 2-channel realization of the multichannel 

model. 

The qq Hamiltonian is embedded in (|T^, but in this matrix formulation we deduce the 
properties of the resonance R due to both its qq interactions and to its AB couplings. Thus, 
finite lifetime effects are included, a priori, in the properties of the qq resonances. Solutions of 
the eigenvalue equation (|^) can be recovered from (|T^) by allowing the off-diagonal potentials 
Vab<^r to become vanishingly small, effectively decoupling R from AB. We demonstrate this 
below. 

All we need do is solve equation (|l^ for the wavefunctions Uas('"ab) and UR^-{rqq)\ This 
will provide us all the information discernable about the AB ^ R system. Describing the 
solutions of equation (|^) is left to Appendicies A, B, and C. 

From the solutions we extract the energy dependent AB phase shift, 6'^'^^{Eab), (see 
equation ( p33| )) due to the presence of the resonance R, and from 6'^'^^{Eab), the mass and 
width of R. We define the width as the energy difference between the energy at 6'^'^^ = 135° 
and 6'^'^^ = 45°, and the resonance mass as the energy at which 6'^'^^ = 90°. 
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A standard description of the phase shift for the process AB —>■ R ^ AB is the non- 
relativistic Breit-Wigner given by 

where and are, respectively, the Breit-Wigner mass and width of the resonance. 
Choosing the quark exchange potential Vab^ab = ^ in equation (H), and therefore ([l^), leaves 
s-channel R production as the only interaction mechanism available to A and B. Therefore, 
the resulting two-channel phase shift, given by equation ( |C33|) in Appendix |C|, should closely 



approximate the Breit-Wigner phase shift if the multichannel model is to be viable. 

In FIG. 6. we plot the two-channel and the Breit-Wigner phase shifts for s-wave Kti 
elastic scattering through the Kq scalar meson. Their excellent agreement establishes this 
essential and non-trivial connection between these two approaches. 

In solving this toy problem we picked a Vab^r coupling strength g that gave the Kq a 



width of only 55 MeV, rather than the physical value of 287 ± 10 MeV |T2[. As we shall see 
below, phase shifts found from scattering through wide resonances deviate from the Breit- 
Wigner lineshape. Also note that the Kq phase shift now passes through 90° at 1.412 GeV, 
down from the eigenenergy value of 1.430 GeV quoted earlier. This mass shift is due entirely 
to the coupling between the Kq and its production and decay Ktc channel, and is a general 
feature of multichannel systems. 

Another required, and non-trivial, connection between this model and the naive quark 
model is that they share spectroscopies in the limits of narrow resonances and no quark 
exchange. In FIG. 7. we plot the Kn phase shifts for different values of g. At the smallest 
value of g the Kq resonant state has a mass of 1.429 GeV and a width of 6 MeV, in excellent 
agreement with the zero width eigenvalue of 1.430 GeV, demonstrating that this correspon- 
dence requirement is satisfied. (Were we to make the Kq still narrower, the scattering energy 
would get still closer to the eigenenergy.) 

Also shown in FIG. 7. are phase shifts at larger values of g. As the lifetime of the 
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resonant state R decreases its' mass also decreases and the Breit-Wigner lineshape is lost. 
Above a critical value of in this case at about g = 1.40 GeV, the Kq mass falls below 
the Kir threshold, thereby making the Kq a linear combination of qq and bound s-wave Ktt 
components, and stable against strong decay. Again, this is entirely a multichannel effect. 

In FIG. 8. we plot sin^(5|f^) for Kit Ko s-wave elastic scattering, with g chosen 
to generate a Kq width of 280 MeV, within the errors of its experimental value. We now 
find the mass of the ground state Kq is 1.355 GeV, a full 75 MeV below the eigenenergy 
of the state in the zero width approximation. Obviously the Hamiltonian parameters need 
to be modified for the multichannel model. The eigenvalue approach can be thought of as 
generating the masses of "bare" qq states; while the multichannel approach "dresses" the 
bare states with connected two-meson channels, making them lighter. The data used to fit 
the bare Hamiltonian parameters has necessarily been the dressed qq spectra, that is all we 
have access to experimentally. This mismatch needs to be recitified; the dressed Hamiltonian 
parameters must be tuned to the dressed data. 

In FIG. 8. we take the Kit scattering energy all the way up to 3.0 GeV, and see the first 
two radial excitations of the Kq at 2.020 and 2.553 GeV, again showing the close connection 
between scattering and spectroscopy. The masses of the first radial and the second radial 
states are, respectively, only 27 and 12 MeV below their eigenenergies. 

It is interesting to explore the behaviour of the Kq resonance produced in s-wave Kn 
scattering. We made a movie of the energy evolution of the scattering and qq wavefunctions 



but were unable to draw any conclusions from it |T3[. In FIG. 9. we plot the energy dependent 
probability of finding some qq state with Kq quantum numbers at the core of the scattering 
system. This probability is given by 



Pko{Ekw) = / \ukoUko\ dr , (16) 



where the ukq normalization is fixed by setting the ampltitude of the external Kn wavefunc- 
tion equal to unity (see Appendix 0). We observe the surprising result that this probability 
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never vanishes. In a beam-beam experiment, at phase shifts of integer multiples of vr radi- 
ans, we find vanishing cross sections and therefore expect that no s-channel qq resonance is 
produced. We have, however, constructed a model in which the incident state is a shell of 
inward falling Kt[ probability amplitude, and not two colliding beams. In this picture it is 
easy to visualize an outward moving shell being phase shifted from an inward moving shell 
by integer multiples of tt radians due to transistions to central qq states. 

We can fix the J of the Kj state by fixing the vertex coupling and the external states, 
thus breaking the Kq, Ki, K2 degeneracy which is present in our eigenstate Hamiltonian. In 
FIG. 10. we show the line shapes for the three processes {Ktt)s ^ Kq, {K*7[)s ^ Ki, and 
{K7c)(i ^ K2 scattering, where the subscript on the two-meson states specifies their relative 
orbital quantum number. We see that the Kj now appear at different masses. The K2-KQ 
mass splitting of 22 MeV arises from their d-wave and s-wave Ktt couplings respectively. 
The Kq-Ki mass difference of 17 MeV arises from Ktt versus K*t[ couplings. These splittings 
are of the same magnitude as the spin-orbit splittings in the eigenvalue quark model [|I|. 

The above discussion also suggests that, when two resonances with the same quantum 
numbers are found at nearby masses in different channels, the correct interpretation is likely 
that they reveal the same underlying state, and that the mass differences are due entirely 
due to different couplings. Testing this conjecture in all resonance channels is an essential 
step in proving the veracity of the multichannel approach. We give an example of this effect 
in the T^tt, Kt], Kt]', Kq system below. 

It is important to emphasize that we are currently only able to demonstrate the qualitative 
behavior of the multichannel solutions. Specific result will require the implementation of 
more accurate models for the AB ^ R verticies, the inclusion of all two-meson channels in 
each sector (which will require the multichannel model discussed below), and a re-tuning of 
the Hamiltonian parameters to some small subset of the low-energy data. While quantitative 
predictions are not yet available the essential points are clear, obvious, and intriguing; s- 
channel meson resonance properties are significantly modified by their two-meson couplings. 
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This multichannel model is an excellent vehicle to begin this exploration. 



V. THE ENERGY DEPENDENCE OF THE BREIT WIGNER PARAMETERS 

Equation ( |(J3ij| ) gives us the energy dependent phase shift in the two-channel model. We 
can use the resulting curve, S^^^lE), to find the energy dependence of the Breit-Wigner mass 
and width parameters of the resonance R. 

Suppose we find 6^'''"{Ei) and 5'^''^{E2) at two neighboring energies Ei and E2. Equating 
^2ch^ from equation ( |C33|) , and 5^*^, from equation (|15D, at these two energies leads to 

If If 

Ci = tan((5^^V^i)) = -r^^ and C2 = iMS'^^E^)) = tt^-^ (17) 

from which it follows that 

E2 C2 - El Ci 



M 



C2-C1 
and 

r« = 2 (M« - El) Ci . (19) 

Using these energy dependent parameters in the Breit-Wigner phase shift equation would 
reproduce the two-channel phase shift for arbitrarily wide resonances. The experimental 
implications are simply that the amplitudes used to represent intermediate isobars must allow 
for variable isobar masses and widths in different regions of the Dalitz Plot. In FIG. 11. and 
FIG. 12. we plot y\.ii{E) and Tii{E) for the two-channel {Kt{)s ^ system as an example 
of their behaviour. Again, specific predictions are required on a channel by channel basis. 

The idea of a variable resonance width has been discussed from the point of view of final 
state interactions in, for example, the reaction t —>■ rj-n-K /r]KK ||T^. Here we see that the 
changing mass and width of the resonance states need not depend on final state interactions 
per se, they arise simply from resonance couplings to final states. 
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VI. QUARK EXCHANGE EFFECTS 



We now consider the possibility of mediating to two-meson scattering by t-channel quark 
exchange effects. 

We omit t-channel meson exchange for three reasons . Meson exchange is topologically 
equivalent to s-channel resonance production and therefore already, at some level, included, 
it is higer order in Ostrong than t-channel quark exchange, and the range for meson exchange 
is of order h/m^^^^,^, which is less than 1 fm for all mesons except the vr. At these separations 
two hadrons, each 1 fm in diameter, have wavefunctions which overlap significantly, and 
quark exchange seems highly probable. 

Much effort has gone into explaining hadron scattering in terms of meson exchange mod- 
els, however, and it would be foolhardy to dismiss this possibility outright, especially at 
larger momentum transfers. In fact, by time-reversing A and C in FIG. 4. we can see that 
t-channel quark exchange can lead to qq exchange diagrams, which may mock-up meson ex- 
change, @J^. Fortunately the multichannel model again offers an excellent vehicle to study 
the inclusion of these terms; they will add to the non-resonant potentials Vab^cd, which are 
generalizations of Vab,-,ab in equation @. 

We were originally led to realize the importance of quark exchange in a variational cal- 
culation of the qqqq ground state wavefunction, ipQqqqq, in the J-^ = 0+ sector @JT^. 'ipoqggq 
required symmetry terms, embodying quark exchange, and the short-range hyperfine in- 
teraction, to represent interacting mesons. Without either feature ipoqqgg represented non- 
interacting mesons. 

The qqqq Hamiltonian used to obtain ipogqgq is just a generalization of the qq Hamiltonian 
of equation (|]), but has added to it a weak color independent quadratic potential, 

= E 1^4 ' (20) 

between all constituent pairs ij that prevents repulsive and non-binding two-meson systems 
from drifting beyond their interaction range. We interpreted the wavefunction 
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4^AB{rAB) = {-ipA S{rA -rB+ rAB)\'ipOqqqq) (21) 

as the amplitude for finding meson A at with wavefunction ip^, and meson B at with 
wavefunction ipg, in ipoqqqq and separated by distance r^B- Inverting the radial Schrodinger 
equation for ip^B allowed us to extract an effective central AB interaction potential Vab^ab'- 

r d2 



2 Ha 



^ABirA 



^-;(r,,) -A'-erj, (22) 



where E^g is the binding energy of the interacting mesons. These led to Gaussian-like poten- 
tials which can be used in a Schrodinger equation, either directly or as "equivalent" square 
well potentials, to find, for example, the phase shifts for 1=3/2 Ktt and 1=2 tttt scattering 
p. These two processes have no s-channel resonances to mediate their interactions, so they 
must go through either t-channel quark exchange, t-channel meson exchange, or some com- 
bination of the two. The phase shifts resulting from the t-channel quark exchange potentials 
compare favorably with experiment |T6|jr7[] . SU(3) relations also allow quark exchange po- 
tentials to be deduced for off-diagonal processes (such as tttt —>■ rjif) from the SU(3) diagonal 
potentials HJT5[ . 



Based on the successes of this technique Barnes etal. [|l^,|T^ developed a more widely 
applicable Born approximation technique for extracting potentials based on one gluon ex- 
change followed by quark (or antiquark) exchange, with interactions normally mediated by 
the short-range hyperfine interaction. For the particular case of meson-meson scattering 



Swanson [19| has shown that the effective intermeson potentials can be modelled accurately 



by the multi- Gaussian form 

Vab^Utab) = E a.e-^(^--/''')' , (23) 

i 

where the and hi depend on the particular reaction AB CD. Agreement between 
these two theoretical approaches and experiment is encouraging, in those sectors where 
comparisons have been made [^^-|TB[. 
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We again choose Kti scattering through the Kq scalar resonance to demonstrate the 
effect of these exchange potentials. First note that the Ktc ^ Kq system is also coupled 
to the Kt] and Krj' systems, so this is a four-channel problem containing three two-meson 
channels and one qq channel. Techniques for solving these multichannel equations are given 
in Appendix 0. 

In FIG. 13. we plot LASS data for 1 = 1/2 Kti scattering ||20| with (a) K-k s-wave 
phase shifts due to two-channel Kti ^ Kq scattering with no quark exchange potentials, 
(b) Ktt ^ Kt] ^ Krj' ^ Kq scattering with only s-channel couplings to the Kq, (c) K-k ^ 
Kr] ^ Krj' scattering with only quark exchange potentials, and (d) Kti ^ Krj ^ Krj' ^ Kq 
scattering with quark exchange potentials. 

In this system the ratio of the relative strengths of the annihilation couplings, which 
we can write as gKn^K„ '■ Qkti^Ko '■ 9Kr]'^Ko where these g's include the Clebsch factors in 
equation (p!T|), are VQ/i^'K-k '■ (1 — V^)//^Xr? : + V^) / l^K-q' 0- The t-channel potentials 
are taken from the variational calculation reviewed above. 

The two-channel phase (a) undershoots the low energy data but fits it in the mid-ranges. 
The phase shift due to only s-channel resonance coupling to the three two-meson states (b) 
undershoots the data below 1.2 GeV and overshoots it above that energy. The t-channel 
quark exchange phase shift (c) looks like an "effective range polynomial background" term, 
suggesting that the need to add a non-resonant background term when carrying out exper- 
imental amplitude analyses is actually a reflection of underlying t-channel quark exchange 
dynamics. The full phase shift (d), including both s- and t-channel processes, is the only 
curve that fits the low-energy data, though above 1.2 GeV it overshoots the data. 

Recall that we are using "bare" Hamiltonian parameters, and a very naive qq annihilation 
and creation vertex model with the range chosen arbitrarily at r,, = 3.0 GeV^^, so the 
deviation from the data is not a statement that the model itself is inadequate, but rather a 
statement that more work is needed to make these qualitative results quantitative. 

As an example of resonance mass shifts seen when an initial two-meson state scatters 
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through a multichannel system into different finals states, we consider 




Kit 



Kr] 



(24) 



In FIG. 14. we show sin^(5' 



;4c/i 

Kit — ^Ktv 



) and sin^(5^ 



■Ach 

K-n^Kri 



) via these four-channel intermediate 



states, using both s-channel Kq couplings and t-channel quark exchange. The apparent Kq 
mass differs by about 45 MeV between these two reactions. 



The multichannel quark model approach treats meson spectroscopy and meson-meson 
scattering on an equal footing. It reproduces the Breit-Wigner and meson spectroscopy 
results for s-channel resonances in the limits of narrow resonance widths and no quark 
exchange. The effects of strong resonance couplings to open decay channels and quark 
exchange effects lead to significant new physics away from these limits. 

This approach must be applied to all of meson spectroscopy and two-body scattering to 
fully test it as a model, and to extract from it the maximum amount of information and 
understanding. 

Of primary importance is the incorporation of realistic AB ^ R vertex potentials, which 
might even be energy dependent functions that would be easily incorporated into the mul- 
tichannel code 0. Each two-meson sector must be considered, and all two-meson and 
resonance states must be included. Following this, the Hamiltonian parameters must be 
refit to some low energy data. At this point the model can be applied to the standard qq 
meson spectrum, where inconsistences abound, as well as to studies of glueball signatures, 
hybrid signatures, and different annihilation models. 

This version of the multichannel model is limited by its non-relativistic kinematics, its 
treatment of external mesons as stable, and its inability to model two-step processes such 
as 



VII. CONCLUSIONS 
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AB 



CDE (25) 



{Rde) C 

y {Rec) D j 

where A, B, C, D, and E are stable mesons, and Rcd, Rde, and R^c are wide and interfering 
intermediate resonances. Developing the model along these lines seems warrented. 

The numbers of open avenues of exploration are numerous. We have pointed out some 
of these along the way, and certainly others remain to be discovered. 

To facilitate the development of the model the fortran source code written to solve the 
multichannel equations is available over the World Wide Web [0]. This code also includes 
some features that were not mentioned in this paper. For example: 

• It already contains a search routine (called fitham.f) to refit the Hamiltonian param- 
eters. This routine calls the qq eigenfunction program (called meson. f) to find the 
eigenenergies of the hadronically stable hadrons, and the scattering program (called 
ccp.f, for coupled channel program) to find the mass and width of hadronically decay- 
ing resonances. With access to this information it searches the Hamiltonian parameter 
space for the best fit to the experimental data. Constraints such as the value of the 
linear confinement potential, or the expected quark masses, can easily be implemented 
as part of the fit. 

• Both the meson eigenvalue program and the scattering program have options of print- 
ing wavefunctions, so the energy evolution of the multichannel wavefunctions can be 
explored |T^. In particular, the contents of the scalar /o(975) and ao(983) KK molecule 



states can be studied once the dressed Hamiltonian parameters are found. This will 
inevitably reveal that, for example, the /o(975) is a mixture of vrvr, KK, rji], t]-)]', 
rj'i]', ^J^{uu + dd), and ss, and that the wavefunctions of each component will vary 
dramatically in shape and relative strength even within the narrow width of the /q. 



18 



• It contains a production model for studying processes such a.s J/ip (/^vrvr and J/ip 
(pKK which hkely proceed through the intermediate 0/q state. The outgoing tttt 
or KK pairs resuh from the creation of a central ss scalar state, which subsequently 
bubbles out through the multichannel scattering potentials into the final state. For 
J/ip ^ UJ7T7T and J/ip uoKK one simply creates a central ^J^luu + dd) scalar state, 
and so on. These processes contain no incoming two-meson states. Solutions are 
found by solving the homogeneous and inhomogeneous multichannel equations, and 
then forming the linear combination of these solutions which have no incoming two- 
meson components. The inhomogeneous multichannel equation is formed by adding 
a 6{0) source term in one of the qq channels to the right hand side of equation (Dl). 
This has the effect of altering the boundary conditions on u(0) and u'(0), but leaving 
the rest of the equation unchanged. The algebra is left as an exercise to the reader. 

• There is a sector for studying baryon-baryon scattering, in particular looking for the 
H-dibaryon and the deuteron as a multichannel effect. 

The multichannel quark model is clearly in its infancy. It is full of possibilities, and holds 
the promise of providing many insights into hadron spectroscopy and scattering dynamics. 
It may prove to be an invaluable, perhaps even an essential, tool for resolving the meson 
spectrum, understanding the qq annihilation verticies, and positively identifying the long- 
sought glueball and hybrid states. In the process it will be teaching us new ideas about 
hadron physics. 
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APPENDIX A: THE NOUMEROV TECHNIQUE 

We introduce the Noumerov technique ^ used to numerically iterate the Schrodinger 
equation 

u"{r) + V{r)uir) = Euir), (Al) 

2jU 



with 



dr dr 

By defining a modified "potential" function W as 



h 

we can rewrite flAll) as 



W{r) ^-'^{E- Vir)) (A3) 



u"{r) = W{r)u{r) , (A4) 

We make r discrete, label the points r^, and set e = rj+i — rj. Taylor expanding M(rj+i) 
and M(rj_i) leads to 

«(r.+i) = «(r,) + en'(r,) + ^,e'u"{ri) + ^,e'u"'{n) + j,e^u""{n) + ... (A5) 
and then to 

^z(r,+i) + n(r,_i) = 2M(r,) + ^"(r,) + ie^M""(r,) + 0{e') , (A6) 

which depends only on even derivatives of u. 
The Noumerov function z{r) is defined as 

z{ri)=u{ri)-^e^u"{r,). (A7) 

By ( |A4D we can rewrite this as 
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z(n)^{i-^,e'w(n))u(n). (A8) 

Note the O(e^) term drops out of the following Taylor expansion 

z{ri+i) + z{ri-i) = {u{ri+i) + - ^e^ + 

= 2 «(r,) + u"{ri) + «""(r,) + 0{e') 

= 2 (u{n) - ^,e^u"{n)) + e^u"{n) + 0(e6) 
^2z{ri) + e'u"{ri) + 0{e') 

^2z{n) + e^W{n)u{ri) + 0{e^) (A9) 
so that, including terms to O(e^), 

z{ri+i) = 2 zin) - z{ri_i) + Win) uin) (AlO) 

and, finally, 

«(r,+i) = (l - ^e^ W^(r,+i)) ^(r,+i) . (All) 

Thus, the Noumerov prescription leads to a value for u{ri+i), given ^(rj), z{ri-i), and 
the VF(rj), which is accurate to O(e^), and which requires a knowledge of the potential V{r) 
only at the points r,. 

To start this numerical proceedure we simply choose u{0) and ■u(e) according to the 
appropriate boundary conditions and then iterate with step size e. As ip{r) oc u{r)/r is the 
solution of the 3-D Schrodinger equation, and finite at r = 0, we require it(0) = 0. The 
choice of u{e) is arbitrary, it merely sets the overall magnitude of the wavefunction. We use 
u{e) — e^"*"^, where £ is the orbital angular momentum, to get the "correct" behaviour of 
^'(0). 
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APPENDIX B: FINDING THE PHYSICAL SOLUTIONS 



In the last Appendix we described the Noumerov technique for solving the Schrodinger 
equation. The physically allowed solutions must obey the additional constraint that, as 
r — > cxD, we must have u(r) — * 0. In this Appendix we describe a short-cut for finding these 
physical solutions with two different potential functions V{r). 

Consider a quantum particle R of mass m bound in a three dimensional central potential 
V{r). The radial equation that must be satisfied is equation (||) 



r t2 



r i{e + 1) 



V{r] 



ui{r) = Eui{r), (Bl) 



2m 

with the boundary conditions ^^(O) = u^(oo) = 0. Here we will only be concerned with the 
large r behaviour of the wavefunction, so the centripital barrier term can be ignored. 
We can, therefore, choose £ = and drop the subscript i. on u^{r) without restricting our 
conclusions. 

Consider first the simple, but important, case of a central square well potential defined 

by 

f Vint if r < rc 
V{r) = (B2) 

I Vext if r > Tc 

where, for the moment, we restrict E to satisfy Vint < E < Vext- 
For r < Tc, ■u(r) = Uintij) obeys 

2^72 

u'Lti'r) = -^iE- Vnt) Uint{r) 

= -kfntUint{r) (B3) 

where the momentum 



/ 2 Tfl I 

',,nt = y - y^nt) = ^J-W{r) (B4) 



is, by construction, real. This has the general solution 
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Uint{r) = Dint sm{kintr) + Gm cos(fcj„tr) . (B5) 

The reason for choosing the unusual nomenclemature, Dmt and Gmt, for the unknown coeffi- 
cients will become apparent below. The boundary condition Uint{0) = automatically leads 

to Gint = 0. 

For r > Tc, u{r) = Uextij) obeys 

= +i^lxtUext{r) (B6) 

where 



/ 2 m 

l^ext = J ~ ^'^^^ 

is, again by construction, real. This has the general solution 

= D,,td{r) + Gextgir) . (B8) 

Analytically, one solves this problem by setting Uintijc) = Uextijc) and u^^^irc) = u'^^tijc) 
to find the unknown co-efficients Dint-, D^xt, and Gext- The physically allowed solutions are 
those for which Gext = 0. It turns out that the values of the D and G depend on the energy 
E through hint and Kext, so the determination of the eigenvalues and eigenfunctions simply 
becomes a matter of finding those E^s for which Gg^t = 0. 



^ We introduce the functions d{r) and g{r) for future convenience. We will continue to express 
the large r wavefunctions of bound systems as linear combinations of some exponentially decaying 
function d{r) plus some exponentially growing function g{r). We will re-define the functional form 
of d{r) and g{r) as we go from the bound square well problem to the linear confining potential, 
but they will play identical roles in all bound cases. For unbound systems d and g will represent, 
respectively, the sin and cos functions. 

24 



Starting with a trial energy of = (Vint + ^E), where AE <^ {Vext — yint), and iterating 
to find u{r) numerically as discussed in Appendix A, one finds that Ge^t is large and has 
the same sign as ^(e). As the trial energy is increased the value of G^xt faUs and eventually 
changes sign. Whenever Gext vanishes we have an eigenenergy and an eigenf unction. Finding 
the zeros of G{E) is a simple matter of interpolation between G{E ) and G{E+), where 
E_ < Egigen < E^, and can easily be done to whatever accuracy the computer allows. We 
would, however, hke to find an efficient algorithm for finding G at each energy. We address 
this now. 

We determine u{r) by stepping out from r = 0, using the numerical techniques of Ap- 
pendix A, until r enters the classically forbidden region. (This region, which begins at r^, 
is the point where the kinetic energy either vanishes or becomes negative.) We record the 
value of the wavef unction at two nearby points > > r^- 

) = Dext d{r^) + GextOira) and 

Uextin) = Dext d{n) + Gext gih) (B9) 

It is then trivial, with u{r) — Uext{i^), to solve for Ggxt 

Gext = ^('^"^ ^('^"^ ~ ^('^"^ (BIO) 
gira) d{n) - g{n) d(rj 

and we can stop the numerical iteration. It is easy to write code to increase the test energy 
by AE until the sign of G changes, then decrease the test energy by some fraction of AE 
until the sign of G changes again, and so on, until one obtains the eigenenergy to the desired 
level of accuracy. In practice we must check to be sure that changing the values of r„ and Vt, 
leaves the eigenenergies fixed. 

When the test energy E exceedes V^xt we no longer have an eigenenergy situation, the 
quantum particle becomes unbound, all energies are allowed, and the external wavefunction 
takes the form 

Uext{r) = Dextsm{kextr) + GextCos{kextr) 
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= D,,td{r)+G,,t9{r) (Bll) 

where 



2m 

ke.t = J ^{E - Ve.M) (B12) 

is real. Notice we have again re-defined the functional forms of d{r) and g{r) without 
changing the role they play in the analysis. Finding u{ra) = Uextij'a) and u{ri^) = Uext{.fh) at 
two different points just outside allows us to find Gext^ via ( |B1(]| ), and D^xt from 

^ ^ u{ra)g{r,) -u{rb)g{r,) 
<i(rj g{rC) - d{r,) g{r,) 

We can now find the phase shift 6, due to the scattering of an unbound particle from the 
3-D square well potential, by noting that 

) = Dext sm{kextr) + G^xt cos{Kxtr) 
= Nsm{kextr + 6) 

= N [sin((5) cos{kextr) + cos{6) sm{kextr)] (B14) 
(with a normalization factor) from which 

5 = atanf^V (B15) 



Thus we use the same algorithm for finding the wavefunction of a bound or scattering 
state; we must simply define the d{r) and g{r) by ( pSD or ( pil|) respectively! For bound 
states we must tune the energy so that Gext = 0; whereas for scattering states all energies 
are permitted and we use ( pi5|) to find S{E). 



In practice this proceedure is found to work beautifully, and the fortran code developed 
to carry out this algorithm runs quickly on a modest unix workstation. 

To solve the more realistic problem of a pure linear confining potential we use the WKB 
approximation P). In the classically forbidden region r > Tc (with Tc = E /h and u{r) = 
Uextij)) we approximate the solution to 
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u"{r) + br u{r) = Eu{r) (B16) 

2m 



by 



|r>rc 



D,,td{r) + Ge^t9{r) (B17) 



where 



<rr) 

'^(E-br) (B18) 



and 



/r 2 771 

V — ^'-^ ~ ^^^^^ 



It appears as though these large r wavefunctions have a different dimension than the 
previous large r wavefunctions (see, for example, equation (pil| )), which would require that 
these expansion coefficients Dext and Gext have different dimensionality that the previous Dext 
and Gext- The discrepancy arises, however, because the sin and cos terms in equation (|B 1 1| ) 
have implicit factors, with dimension 1 / -y/ftn and magnitude unity, that are required to 
normalize the otherwise divergent sin and cos functions to one particle per unit length. 

Notice also that we have again re-defined the functional form of d{r) and g{r) between 
equations ( |B 1 1| ) and ( [Bl?] ) without changing the role they play in our analysis. 

We can again write two nearby values of the wavefunction, as in equations (PP|), and find 
Gext by equation ( |B1CI| ). Knowing Gext makes it easy to find the eigenvalues and eigenfunc- 



tions of equation (|B16| ). 

To solve this problem when the differential equation is given by (^) and (P) we note that, 
in the classically forbidden region, we can make the approximation 
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I2 



+ ^™ (^m) ^m^ + m^ + C + hr (B20) 



2 l^qq Tqq 

which allows us to use the same scheme as in the pure linear case, except that is now 
determined by the condition 



^ ^qqi^^qq ~\~ 1 

2^ ' 



+ K rre)=E, (B21) 



and 



/ 2 m 

K= J-^{E -m^-m^-C -br) (B22) 



while 



/r / 2/71 
W — ^ (-^ — rrig — rrig — C — hr) dx 

2 / 2rr? 

= ^ \/-^(^ - - C - 6r)3 (B23) 

In practice these approximate schemes must be monitored to ensure that the solutions 
do not depend on the numerical parameters e, r^, and introduced. The model and the ex- 
perimental results are also inexact, and the numerical errors are easily rendered insignificant 
by comparison. 

APPENDIX C: SOLVING THE TWO CHANNEL MODEL 

In this Appendix we discuss techniques for solving equation ( p!4D for the specific case of 
AB ^ R ^ AB. Here A, B, and R are ordinary qq resonances, but A and B are assumed 
to be stable against strong decay while R obviously is not, it can have a lifetime of 0(10^^'^) 
seconds. 

Our first step is eliminating one of the variables r^B and Tqq. We take the SU(3) limit of 
Tqq = 2 Tab- This substitutiou is easily implemented by writing equation (|l^ in the form, 
with r = r4B, 

28 







= E 















(CI) 



VAs^nir) H,j2r) 

We further simplify this equation by replacing the qq wavefunction Ua^-^{2r) with Ua{r) = 



(C2) 



We can write the above equation in a more compact form if we first define the Hamiltonian 
matrix 



HAB{r) VAB^nir) 




UAsir) 




UAsir) 








= E 




VAB^nir) H,j2r) _ 











n(r) 



and the vector radial wavefunction 



u r 



(C3) 



Mi(r) 




UAsir) 




UAB{r) 


M2(r) 




_ ««,,(2r) _ 







(C4) 



Explicitly, 



+ mg + mg + Vn (2r) 



(C5) 



and, from equation (^, 
(2r) 



3 etc 0" 
- 2r — 

4 2r ^ 7r3/2 



3 
4 



F F- 



3 rrigmg vr^/^ g g g g 



(C6) 



With all of this equation (|Cl| ) becomes 

Ti. (r) u(r) = Eu{r) . 



(C7) 



From equations (^, (|^), and (|TT]), we see that each component of H consists of a function 
of r plus, in the diagonal entries, a derivative operator. We define a new matrix W, whose 
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entries are proportional to these functions of r, and whose two indicies refer, respectively, to 
the final and initial channels, by 



W(r) 



Wu{r) Wu{r) 
W2i(r) W22(r) 



(a 



with channel 1 being the AB system and channel 2 the qq state R. Explicitly, 



(r) + mA + mB - E 



Wo^(r) 



h 

2 fJ'qg 



f h 1 + 1) ^ +rn, + m,-E 



2'V^B^Jr) 



' h ' ' IJ'AB 

The factors of 2^ in (|C9D effect the change from Tqq to 2r. 
We can now write equation ( pl|) as 

u"(r) =W(r)u(r). 



(C9) 



(CIO) 



( |C10|) is equivalent in every way to equation ([T^), but it lends itself to a matrix formuation 
of the Noumerov technique and is easily adapted to generate numerical solutions. 
Analagous to Appendix 0, we define the vector Noumerov function as 



z(r) = u(r) — j^e^ u"(r) 



(Cll) 



so, by dClOD 



z(r) = (j-ie^ W(r))u(r) 



(C12) 



where X is the unit matrix. 

When we make these functions depend on the discrete variable and Taylor expand 
u(rj) we find, in analogy with equation ( [A6|) , 
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u(r,+i) + u(r,_i) = 2 u(r,) + u"(r,) + ^e^ u""(r,) + 0(e'= 
and, in analogy with equation ([A10| ), 

z(r,+i) = 2 z(r,) - z(r,_i) + >V(r,) u(r,) + O(e^) . 



(C13) 



(CM) 



To find u(rj+i), which is the wavefunction we are interested in, we simply invert equa- 
tion ( ICT2D : 

u(r,) = (j -ie^ n;(r,))"' z(r,) (C15) 

Therefore, knowing u(rj), z(rj), z(rj_i), and the W(ri), allows us to calulate u(rj+i): z.e., 
we have an algorithm for finding u(r), given u(0), u(e), and VV(r), which is exact to O(e^). 

Having found the mathematical solutions to equation (|C1D we must now, as in Appendix 
0, find the physically allowed solutions, namely those which have Uqq{rqq) vanish as Tqq — > 00. 

As in the case of the single channel Schrodinger equation, the iterated numerical solution 
will have, generally, both exponentially growing and decaying components in the classically 
forbidden region. In the single channel case we were able to adjust the test energy of the 
bound state to find the eigenenergy and the physical wavef unctions. In the two-channel 
problem, however, the energy is set by the total energy which A and B bring to the interac- 
tion, and has any value greater than + m^. We must, therefore, find a new criterion for 
defining the physical wavefunctions. 

The two-channel large r wavefunction is 

Di sin(A;ir) + Gi cos(fcir) 



D2 e-y^^^y,/^ + G2 e+y^^^y^^^ 



(C16) 



with 



and 



y2[r) 



-{E - itla 



■—^{E — nig — nig — C ~ br) 



2 

36 



2/i2 



{E — rriq — nig — C — bry 



(C17) 
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and Tc defined as the point at which the qq kinetic energy vanishes, corresponding to the 
boundary between the classically allowed and forbidden regions. The upper component 
describes two free mesons with kinetic energy {E — tua — rrig) and the lower component, 
written in analogy with equation ( |B17| ), describes a qq state bound by a linear confining 
potential. The unknown coefficients Di, D2, Gi, and G2 are determined using the techniques 
discussed in Appendix B. 



It is obvious that the solution we want has G2 = in (|C16| ). Suppose we generate two 
different solutions of equation ( |C10|) according to two different sets of initial conditions: 



ui(0) 








and Ui(e) 



e{^i+i) 




(C18) 



and 



U2(0) 








and U2(e) 





.(^2 + 1) 



(C19) 



where ij is the orbital angular momentum of the j*'' channel. These will generate the 
aymptotic solutions 



Dii sin(A;ir) + Gu cos{kir) 
D21 e-^^M/ + G21 e+^^M/y^ 



(C20) 



and 



U2[r) \r>T 



(C21) 



D12 sin(A;ir) + G12 cos{kir) 
D22 e-y^^^y^^) + G22 e+^^('^V\A^ 

where we have added an index to the D and G relative to equation ( |C16|) to accomodate our 
having two solutions. By defining a matrix of wavefunctions as 



U (r) = [ui(r) U2(r)] , 



(C22) 



( |C20D and (|C2lD can be written 
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U (r) 



k>rc 



sin(/cir) 

e-V'^''^ J K2{r) 



+ 



For future convenience we define the matricies 



cos(A;ir) 

e+y'^'-y J K2{r) 



D21 D22 

Gn G12 
G21 G22 



(C23) 





Dn 


D12 




Gil G12 


V = 












D21 


D22 




G21 G22 



and C = Q ^ 



Cll C12 
G21 G22 



(C24) 



where we know that C exists because Ui(r) and U2(r) are hnearly independent solutions. 
By noting that 



{n {r)-IE)U=0 
we see that UC is also a solution of the two-channel Schrodinger equation: 

{n {r)-I E) UC=0. 
Using equation ( |C23| ) we can write 





(C25) 



(C26) 



U (r) C 



\r>rc 



+ 



sm{kir) 

e-y^'^''^ ^ K2ir) 

cos(fcir) 

6+2^2 M/y'/talr) 



1^11 D12 
D21 D22 



Gn G12 
G21 G22 



(C27) 



We note that the first column of ( P27|) is exactly the solution that we are seeking, it has no 
exponentially growing component! 

It is useful to simplify this equation by introducing yet another matrix: 



Fll F12 




Dn D12 




Gn G12 


F21 F22 




D21 D22 




G21 G22 



--VC 



(C28) 



where 
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kj ■ 



(C29) 



k=l 



The physical solution Up(r) is the first column of the solution matrix U{r)C: 



Up r 



\r>rc 



\U{r) 



r>rc 



Cli =N 



Fii sin(A;ir) + cos(A;ir) 

where is a normalization factor. This solves equation ([T^)! 

If we normalize the AB wavefunction to unit amplitude (per unit length) then 



(C30) 



(C31) 



where 



^21 



FK + 1 



(C32) 



and the 2-channel phase shift induced as AB forms R and R then decays to AB is 

1 



5 



2ch 



at an 



'11 



(C33) 



By solving for 5'^'^^{E) at a range of energies beginning at AB threshold we can compare 
the predictions of the coupled channel equations directly with the canonical Breit-Wigner 
phase shift for AB ^ R AB. This comparison is made in FIG. 6. and discussed in 



Section IV 



Since 



W(0) 








and 



Uit] 



e(^i+i) 



(C34) 



and 



Up(r) = [W(r)C]i, (C35) 

we can generate these physical solutions numerically by starting with the intitial conditions 

34 



Up(0) = [W(0)C] 























- 1 





(C36) 



and 
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e(^^+i) 




C21 
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t-i) 



(C37) 



We can numerically integrate the lower component squared of equation ( |C3 1| ) to obtain 
the probability of finding the AB ^ R system in R. A plot of this probability as a function 
of the scattering energy E is shown in FIG. 9. Note that the probability of R being present 
in the central region never vanishes, even when the phase shift of AB-in to AB-out equals 
an integer multiple of vr. This is because AB-in is a shell of inward-moving AB probability 
amplitudes rather than two columnated beams. 



APPENDIX D: SOLVING THE MULTICHANNEL MODEL 

At this point we turn to a more general multichannel problem that contains essentially 
all the special cases we must consider; the two-channel equation (CI) is generalized to the 
four-channel equation for AB ^ CD ^ Q ^ R: 





Vab^cd (''") 


VAB^air) 


VAB^nir) 




UAB^r) 




UAB{r) 




Hco{r) 








Ucoir) 


= E 










M,^Bir) 








Mr) 




VcD^R{r) 


M,^B{r) 


Hnj2r) _ 











(Dl) 



Here AB and CD represent two different meson-meson final states, with m^+mg < mc+rrioj^ 
and Q and R represent two meson resonances with tttq < m^. 



^ As an example, to study the mass splitting caused by different decay channel couplings, and 
the D/S amplitude ratio, for the reaction (lott)^ <^ (w/r)^^ •s-^ 6i(1235) we would construct a three- 
channel equation with + itib = fric + m^,. 
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In (Dl) those elements already introduced are Hj,b and Hen, given by equation (|^), 
or its analogue, Hq^^ and H^^^, given by equations (|^) and (|^), or their analogues, and 
Kis«Q5 Vab*-,r, Vcd^q, and Vcd^r, given by equation (p^l]), or its analogues. There are two 
new elements, Vab^cd is the quark exchange potential reviewed in Section |VI|, and Vq^r, 
represents possible direct meson-meson mixing such as 77 — 77' or — mixing, which we 
set = here. Allowing this entry to be non-zero is one way to mock-up the presence of a 
glueball state. ^ 

Analogous to the two-channel case (Appendix C), we define a 4x1 vector wavefunction 



u r 



Ui 


[r) 










U2 


[r) 




Ucoir) 
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[r) 
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[r) _ 











(D2) 



and a 4x4 matrix W (r) with entries given by the direct analogues of equations 

We again form a matrix of, in this case, 4 solutions Uj(r), with j = 1,4, and use the 
Numerov technique to find these wavefunction given 4 orthogonal initial conditions: 



W(0) 



Ui(0) U2(0) U3(0) U4(0) 































and 



more realistic way would be to form a separate glueball channel analogous to a qq channel. In 
our non-relativistic formulation this would require modelling the diagonal glueball Hamiltonian in 



terms of "constituent" valance gluons [22 
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e(^^+i) 







U{e) 



(D3) 







e(^3+i) 











e(^4+i) 



where each column forms a separate and linearly independent solution of the four dimensional 



Each of these solutions will have one of two possible large r expansions for the CD 
system, depending on whether the energy E is smaller than or larger than iric + m^,, so that, 
respectively, CD is either bound or free. 

If E' < iric + mj, then CD may exist as a bound state confined to the central region of 
the interaction, but not as a free state. (Recall that we are thinking of A, B, C, and D as 
qq states which are stable against strong decays.) In the large r regions all CD potentials 
have vanished and the CD confinement is due solely to conservation of energy. Therefore, 
when we have found the wavefunctions by iteration from r = to > > using the 
techniques described in Appendix 0, we must represent the large r bound CD wavefunction 
as the sum of an exponentially decaying part and an exponentially growing part, in direct 
analogy with the square well discussion in Appendix p. The AB system remains free and 
the resonances Q and R are bound by linear confining potentials as in Appendix |C[ For the 
j^^ solution these large r wavefunctions are 



form of equation ( CIQ ) and ij is the orbital angular momenta of the j*'' channel. 



Dij sm{kir) + Gij cos{kir) 

D2i e-'^^r + Q e+«2r 




(D4) 




K2 plays the same role that Ke^jt 



played in Appendix |B[ All other quantities are as defined as in Appendix y, and the D and 
C coefficients are found as described in Appendix 0. 



37 



With these four solutions we go through exactly the same proceedure as discussed in 
Appendix y, but now all matricies are 4x4 and T has its 16 components defined by 



(D5) 



k=l 



As in Appendix 0, we find the large r form of the physical solution up to be the first column 
of the solution matrix U (r) C: 



Mr)\r>r. = [U {r)\r>r. C]i =N 



Fii sin(A;ir) + cos(A;ir) 



(D6) 



where is a normalization factor. 

The phase shift of the AB state resulting from AB 
scattering process is 

1 



c-4c/i 



at an 



'11 



(AB ^ CD ^ Q ^ R) ^ AB 



(D7) 



Since we are below the energy for inelastic scattering to CD, all the AB amplitude fiowing 
into the interaction region must also fiow out as AB. We choose, for later convenience, 
to normalize the spherical inward-moving external wavef unction to unit fiux crossing an 
imaginary sphere surrounding the interaction region, which leads to. 

\/4:-kvab (Kem^abm) = V47rf 1 (m*.^ = 1, (D8) 

where Vab = vi = ^2 {E — tjia — rng)/fJ'AB is the "reduced" velocity of AB (channel 1) in. 
This has the effect of implementing conservation of particle number |^3|. With this choice 

sin(/cir + 5^'^'*) 
N2 e-'^^r 



1 



Up r 



\r>rc 



(D9) 
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where, we see from 



that, for j = 2, 3, 4, 



+ 1 



(DIO) 



When E > rric + m^, we are in a very different regime, CD is now a free state and we 
can describe inelastic scattering processes such as AB CD. We use the same proceedure 
as above to find 4 solutions to (Dl), only now we must represent the large r wavefunction 
of the j^^ solution by 



Dij sin(fcir) + Gij cos(fcir) 
D2j sin(fc2r) + C2j cos(A;2r) 



(Dll) 



where the momentum k2 = y {2 ^cd/T^^){E — rric — nio) is real. In this case there are now two 
physically allowed solutions, corresponding to the first two columns of the matrix U (r) C. 



We find, in direct analogy with 



Upi(r)|r.>rc = [U {r)\r>r,C]i 



that these two solutions are 

Fii sin(A;ir) + cos(A;ir) 
F21 mi{k2r) 



F31 e-^^W/\A^ 
F41 e-^^W/^^ 



(D12) 



and 



|r>rt 



c 



(D13) 



F12 sin(fcir) 
F22 sm{k2r) + cos(A;2r) 

Of course, now every linear combination of (|D12| ) and (pi3|) is also a solution of the 
multichannel equation. We are interested in solving for two specific linear combinations, one 
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corresponding to an incident state of pure AB, which models the inelastic process AB AB 
and CD, and the other with pure CD in, for CD —>■ AB and CD. 

For pure AB incident we wish to find the complex coefficients an and a2i so that 

{Fii an + Fi2 8.21) sm{ki r) + an cos(A;i r) 
(F21 an + F22 a2i) sin(fc2 r) + a2i cos(A;2 r) 



(F31 an + F32 a2i)e-^='M/V'«3(r) 
(F41 an + F42 a2i) e-^^W/^^ 

1 e-'^^^ + An e+^^i^ 
e-^'^^r- _^ A21 e+^^^r- 

(F31 an + F32 a2i) e-^^M/ 



(D14) 



(F41 an + F42 a2i) e-^^M/^/€4(r) 
where An and A21 are the complex coefficients of the outward-moving AB and CD systems, 
respectively, and are to be determined from the numerical solution. All effects of the inter- 
mediate resonances Q and R on the scattering states will be realized in the coefficients An 
and A21 of the outward-moving waves. We shall address the issue of normalization when we 
relate these solutions to the ^-matrix below. 

We see that only the upper two components of ( pi4| ) will help us solve for the 4 complex 
unknowns an, a2i. An, and A21. As e^*'^^' = cos{kr) ± ism{kr), equating the coefficients of 
cos(/cir), sin(/cir), cos{k2r), and sin(fc2'") leads to the four complex equations 

Fu an + F12 a2i = -i + iku , 



F21 an + F22 a2i = iA; 
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an = 1 + A 



a2i = A21 . 



11 ) 



and 



(D15) 



By defining 
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an 




An 




, ai = 






F2I F22 




a2i 




A21 



jr(2) 

equations (pi5|) can be written in matrix form 

P''^ ai = -iei + iAi 
ai = ei + Ai . 



and ei 



1 




(D16) 



(D17) 



(Note that, although JF'^' is now a 2x2 matrix, it still depends on U3(r) and U4(r) because 
the sum in equation (P5| ) runs from 1 to 4.) From this one can show that the required 
complex co-efficients are 



ai = 2 [X'^' -i ] [ J(^' + ei 
and the resulting complex amplitudes of the outward-moving waves are 
Ai = 2 [J(^) -I ] [ J(^> + ei - ei 



(D18) 



(D19) 



where X*^' is the 2x2 unit matrix. These completely specify specify the physical wavefunctions 
for AB incident. 

Solving this problem for pure CD incident simply requires that we find the complex 
coefficients ai2 and a22 so that 

UcD;in('")|r>rc = ai2 Upi(r)|r>r, + a22Up2(r)|^>r^ 

(Fn ai2 + Fx2 ^22) sin(A;i r) + ai2 cos(A;i r) 
(F21 ai2 + F22 a22) sin(A;2 r) + a22 cos(A;2 r) 
(F31 ai2 + F32 a22) e-^3(r)/y^ 
(F41 ai2 + F42 a22) e-^*(^'-Vy^ 
e-^'^i^ + A12 e+^*=i^ 



1 e-'^^'' + A22 e+*'=2^ 
(F31 ai2 + F32 a22) e-y--'^'-^ 
(F41 ai2 + F42 a22) e-^^(^V\/«^ 



(D20) 
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where A12 and A22 are the complex coefficients of the outward-moving AB and CD systems, 
respectively. 

The solution to this problem is found by defining 





ai2 




A12 







(D21) 


a2 = 


, A2 = 


, and §2 = 






a22 




A22 




1 





and then by replacing ei in equations (|D18| ) and ( |D19| ) by §2, which leads to solutions for 
a2 and A2: 

a2 = [J'^' -2 J^'^' ] [X<'' + T^'^ J^'^' 62 (D22) 

V47rf2 

and 

A2 = 2 [J(') -i ^('^ ] [ J''> + ^('^ 62 - 62 (D23) 

Since we now know all four wavefunctions for the four possible two-channel scattering 
processes AB ^ CD we can find the two-channel S'-matrix from the parameterization 



Sii 


S12 




r^e^*^! i[l 


_ ^2ji/2 gj(<5i+<52) 
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S22 











To express the parameters rj, 61, and 62 in terms of the coefficients of the outward-moving 
wavefunctions we need to consider wavefunction normalizations and conservation of particle 
flux. 

Notice, from equation (pi4|) , that the external part of the wavefunction consists of a shell 
of inward-moving AB and a shell of outward-moving AB and CD. We can express this as 







{An}e+*'=i^ 


(D25) 




{0} e-^'^^r 




{A2i}e+*'=2^ 





Conservation of particle flux requires that the number of inward-moving and outward- 
moving particles must be equal in any given unit of time. We can realize this by normalizing 
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the inward-moving wave to unit flux, i. e. multiplying through by a factor of 1/ y/Airvi, which 
leads to 



{0} e-*^2r 



{Aii}e+*'=i^ 



1 



V47rt>i 
{An}e+^'=^7v/4^ 



(D26) 



Thus the multichannel quark model takes an inward-moving unit flux normalized AB 
system having an amplitude of 1 into a linear combination of outward-moving unit nor- 



malized AB and CD with complex amplitudes An and A2iyf2/fi respectively. (As CD is 
heavier than AB (by construction), they have lower relative velocity at a given total energy, 
and so the wavefunction describing their relative separation must have greater density per 
unit volume to satisfy conservation of particle flux.) 
For pure inward-moving CD (|D2CI| ) says 
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{Ai2vW^} e+^^^^Vv/i^ 



{A22}e+*'^27v/4^ 



(D27) 



From ( P24D , ( P26D and (|D27|) we find that 
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Sii 
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S21 


S22 





All Al2\Jvilv2 
A2l\Jv2/Vi A22 



(D28) 



so that, with (|D24|) 



-atan ( ^ I . 
2 VMAii}, 



-atan 



5i 



<5i 

?? = ^(Im{Aii})2+(Re{Aii})^ 
which solves (Dl) for two free two-meson states. 



/Re{A2i}^ 
VMA21}, 



and 



(D29) 
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We can also write the cross section, in millibarnes, for an initial state i to scatter into a 



final state / as |23 



(Tfi 



(Re{S^J - 5f,y + (Im{S^J)^J (D30) 



where kf is the momentum of the outward-moving state in inverse fermi and 6fi is 1 when 
f = i and zero otherwise. 

The generalization of this proceedure to arbitrary numbers of external states AiBi, A2B2, 
^3^3, ... and intermediate resonance states is a straightforward exercise left to the reader. 
We will, however, only require only one or two resonance states until we get to the glueball, 
hybrid, or charmed quark sectors: in isospin we have the ^J^{uu + dd) and ss states, while 
in 1=1 and 1/2 we have only one possible qq resonance. To date the largest system solved is 
the seven channel, two resonance, tttt, KK, rji], rji]', ri'i]', /q, /q system describing s-wave 
7T7T scattering |p. 
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FIGURES 
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FIG. 1. Here the Vqq potential has the quantum numbers of the Kj, l^^ = 1, and s^s = 1- 
{Vqg{r) — El) and 25ugg{r) are plotted for two solutions of equation (5), with trial energies 
El = 1.4299235012 GeV and E2 = 1.4299235014 GeV. The actual eigenenergy lies between these 
values. We plot {Vgg{r) — Ei) rather than Vgg{r) so that the zero of n(r) and the zero of available 
kinetic energy both lie on the r axis. We pick the wavefunction normalization of 25 simply so we 
can see the wavefunction and the potential together. 
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r in fm 



FIG. 2. The same as FIG. 1., but for the first radiahy excited Kj state. The trial energies are 
El = 2.047471763 GeV and E2 = 2.047471768 GeV. 
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FIG. 3. A quark line diagram for s-channel meson resonance formation in meson-meson 
scattering. The solid lines are quarks, the dashed lines are anti-quarks, and the shaded dots are qq 
annihilation or creation vertecies. Understanding the properties of the resonance R in terms of the 
underlying qq interactions and the AB R ^ CD couplings is a goal of modern particle physics. 
As the lifetime of R is typically 0(10^^^) seconds it is clear that the properties of R may depend 
significantly on its couplings to AB and CD. 
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FIG. 4. A quark line diagram for one-gluon exchange (the curly line) followed by quark 
exchange in meson-meson scattering. This diagram leads to non-resonant hadronic range van der 
Waals type intermeson potentials which have gaussian shapes and ranges of about 0.5 fm. 
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FIG. 5. The relation between Vab and Vji = Vqg. Here A = qiqa, B = q2qi, and R = qiq^. 
Clearly = 2 Tab when all quarks are of equal mass. 
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FIG. 6. Compare sin'^{d^) (the dots) and sin^((5^^) (the soHd Hne) for s-wave Ktt scattering. 
Here 5^ is the Breit-Wigner phase shift and (5^^ (the sohd hne) is the phase shift found from the 
two-channel quark model, with the quark exchange potential in the quark model equal set to zero. 
The agreement between the two descriptions is excellent. 
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FIG. 7. The variation in the 2-channel Kir s-wave scattering phase shift S'^'^^ as g varies. The 
horizontal dotted hnes mark 45°, 90°, and 135°, and the numbers closest to the phase shifts give 
the value of g in GeV. Note the departure from the Breit-Wigner phase shift description as the 
width increases, and the onset of a Ktt bound state as g grows larger than 1.40 GeV. We define the 
width as the energy difference between the energy at (J^^'' = 135° and 6^"'' = 45°, and the resonance 
energy as the energy at which 6'^'^^ = 90°. 
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FIG. 8. sin^((5^^), from s-wave K-k scattering, is plotted from threshold to Ek^k = 3.0 GeV, 
with g = 0.60 GeV and Tj^^^ = 280 MeV. The peaks at 1.356, 2.020, and 2.553 GeV correspond to 
the ground state, the first radial excitation, and the second radial excitation of the L = 1, 5 = 1, 
J = strangeness ±1 resonance. The narrowing of the excited states depends on the form of the 
qq annihilation/creation verticies and requires further study. Note that heavier states have more 
open channels, and that here we are looking only at the partial width into a single channel. 
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FIG. 9. The probability, /q°° \u*i(^uko\^ dr , that the s-wave Ktt scattering state is in a central 
qq resonance with Kq quantum numbers, plotted as a function of the Kn scattering energy. Note 
that some us is always present in the central region. The location and widths of these peaks are 
directly related to the mass and width of the qq resonances. 
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FIG. 10. Resonance masses of the strangeness ±1 £qq = 1, Sqq = 1, J = 0, 1, 2 triplet as J varies. 
The Kq and K2 are obtained from the s- and d-wave scattering of a Kit system respectively, while 
the Ki is obtained from s-wave scattering of a K*tt pair. The resonances masses are 1.397, 1.370, 
and 1.419 GeV for J = 0, 1, and 2 respectively. The K2 — Kq mass splitting of 22 MeV arises from 
the different vertex couplings to Ktt whereas the Kq — Ki mass difference of 17 MeV arises from 
s-wave couplings to different final states. The underlying Kq, Ki, K2 masses, in this version of the 
naive quark model are degenerate since we have omitted spin-orbit couplings in the Hamiltonian. 
The strength of the vertex couplings have been adjusted so each resonance has the same width in 
order to reduce the number of differences and clarify the discussion. 
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FIG. 11. The dependence of the Breit-Wigner mass as extracted from the two-channel solution 

for Kn ^ Kq scattering. The dots indicate the resonance masses for the ground, first radial, and 
second radial states. There is a third radial excitation just above 3.0 GeV. 
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FIG. 12. The dependence of the Breit-Wigner width as extracted from the two-channel solution 
for Kt: ^ Kq scattering. The dots indicate the width predictions at the resonance energies for the 
ground, first radially excited, and second radially excited states. The location of the dots below 
the peaks indicates that the resonances are narrower on the low side than on the high side. The 
details of these results will change as more realistic models for the annihilation potentials are used. 
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FIG. 13. The Kit phase shift from the four-channel equation, (a) uses only Kir <^ Kq 
couphngs with no quark exchange potentials (so it's just the Breit-Wigner phase shift), (b) uses 
only the resonance couphngs of the three two-meson states to the Kq, (c) uses only quark exchange 
potentials, and (d) uses both quark exchange and resonance couplings. The data is from LASS [20]. 
The two vertical bars denote the onset of the Kr] threshold at 1.045 GeV and the Kr]' threshold 
at 1.445 GeV. See Section VI for a discussion of these solutions. 
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FIG. 14. We plot sin^((5^^_,^^) and sin^((5^^^^^) for the scattering processes Kt: —>■ 
(Kir ^ Kr] ^ Kr]' ^ Kq) Kit and Kir (Kit ^ Kt] Kt]' ^ Kq) Kt] where the 
curves are labelled by their final states. The 45 MeV mass difference in the central Kq resonance 
in these two reactions is a pure multichannel effect. The decrease in sin 

{SK'i!r^Kri) Kr] threshold 

suggests that there is a large Kt] bound state component in the multichannel system below this 
energy. We use both s-channel resonance couplings and t-channel quark exchange to generate these 
curves. 
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